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We develop a method to determine accurately the binding energy of the 
X(3872) from lattice data for the DD interaction. We show that, because 
of the small difference between the neutral and charged components of the 
X(3872), it is necessary to differentiate them in the energy levels of the 
lattice spectrum if one wishes to have a precise determination of the the 
binding energy of the X(3872). The analysis of the data requires the use 
of coupled channels. Depending on the number of levels available and the 
size of the box we determine the precision needed in the lattice energies 
to hnally obtain a desired accuracy in the binding energy. 
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1 Introduction 


The purpose of the present paper is to find a strategy to determine accurately the 
binding energy of the X(3872) in lattice QCD simulations. A precise determination, 
with an energy about 0.2 MeV below the D^D threshold, requires to differentiate 
between the u and d quark masses in order to account for the 7 MeV difference between 
the neutral and charged components of the wave function [1]. The small binding of the 
state with respect to the D^D threshold, much smaller than the difference of masses 
between the D^D and D'^D*~ components, makes this consideration imperative in 
order to get a precise value of the binding energy and unambiguously determine the 
bound state character of the X(3872). In fact, when this is done, energy levels can 
be associated to either D^D or D^D . 


2 The X(3872) in the continuum limit 

In this section we discuss briefly the dynamical generation of the X(3872) in the con¬ 
tinuum limit. All the details are in Refs. BEI. The pseudoscalar - vector interaction 
can be studied through the hidden gauge Lagrangian [3], which contains interaction 
between vectors and with pseudoscalar mesons. The model is based on vector-meson 
exchange, see Fig. [H In fact, for s-wave, when the momenta exchanged in the 



Figure 1: Point-like pseudoscalar - vector interaction. 


propagator of the vector meson exchanged can be neglected against —My, leads to a 
point-like interaction, and is equivalent to using the Lagrangian, 


C 


PPVV 


4f 




( 1 ) 


with = {d^P)P — Pd^P and — ViydfjV'', see [D E]. In Ref. pQ, the 

currents are separated for heavy and light vector-meson-exchange, introducing the 


breaking parameters, 7 = V’ = with mg. = m^ = 800 

MeV, m 3 * = rriH = 2050 MeV and mi* = mj/^ = 3097 MeV. This gives, 7 = 0.14 and 
V’ = 0.07. Because of the smallness of the breaking parameters, the light and heavy 
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sector are almost disconnected, and the transition potential between those is very 
small. Also, for light mesons, / = /tt = 93 MeV, and for heavy ones, f = fn = 165 
MeV, is used. Thus, the amplitude of the process Vi{k)Pi{p) —)■ V 2 {k')P 2 {p'), is given 
by 


Vij{s,t,u) 






u) e.e' 


( 2 ) 


with s — u = {k + k'){p + p'), which must be projected in s-wave [Hll], and i,j refer to 
the particle channels. Working in the charge basis, we have the channels K^ — 

C.C.), - C.C.), - C.C.), - c.c.) and - c.c.), 

and the matrix ^ can be written in this basis as 
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Eq. ([2]) is the input of the Bethe Salpether equation, 

T = {I - VG)~We.e'. (4) 


Here G a diagonal matrix of the two-meson loop function for each channel. Usually 
it is evaluated with dimensional regularization and depends on the parameter a [ 1 ] 
ore one can also evaluate the G function with a cuttoff. The calculation in Ref. |1] is 
redone to get a binding energy more realistic at 0.2 MeV with respect to the channel 
-c.c. m, where the masses of the mesons are taken from the PDG [5]. The free 
parameter, a, is hxed for the light channels, ai = — 0.8 Bill. but the pole position of 
the X(3872) is not sensitive to that, since its mass is far away from these thresholds. 
For the heavy channels, the value an = —1.265 is needed for such binding energy (/i is 
taken equal to 1500 MeV in all channels). In Tabled! a summary of the pole position 
and couplings of the resonance to each channel is given. The Weinberg compositeness 
condition [B] can be generalized for dynamically generated resonances from several 
channels [1] and gives the probability of Ending the i channel in the wave function, 
which are 0.86 for D*^D^ — c.c, 0.124 for D*^D~ — c.c and 0.016 for D~ — c.c. 
However, this is difierent from the wave function at the origin (27r)^A.^(^o)i = QiGi, 
which usually enters the evaluation of observables and are nearly equal [ 1 ]. 


3 Formalism in finite volume 

We follow the formalism used Ref. [7] where the infinite volume amplitude T is re¬ 
placed by the amplitude T in a finite box of size L and G(P°) is replaced by the finite 
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V^o = (3871.6 - iO.OOl) MeV 


Channel 

\g,\ [MeV] 

^(A'-A'+ - c.c) 

53 

- c.c) 

49 

UD'*D- - c.c) 

3638 


3663 

- c.c) 

3395 


Table 1: Couplings of the pole at MeV to the channel i. 


volume loop function denoted with G, given by the discrete sum over eigenstates of 
the box 

Qi 


with 

r(pO +^2(gl) _1_ 

^ 2u;,{ql)u2iql) (P^y - Mql) + U2{ql)r 


( 6 ) 


where cUj = \/rnf + {ql ^ is the energy and the momentum ^is quantized as qi = 
corresponding to the periodic boundary conditions. Here the vector f?, denotes the 
three dimension vector of all integers (Z^). This form produces a degeneracy for the 
set of three integer which has the same modulus. The sum over the momenta is done 
until a qmax, so the three dimension sum over n* in Eq. (|5]) becomes a one dimension 
sum over rrii to an Umax in a symmetric box, Umax = ^2^^. The equivalent formalism 
in hnite volume should also be made independent of qmax and related to a. This is 
done in Ref. [8] with the result 


G 


lim 

Qmax — 





= G^^+ lim 5G. 

Qmax —^00 


(7) 

In Fig.[2]we show that 5G converges as qmax —t 00. In practice, one can take an average 
for different values between qmax =1500-2500 MeV and one sees that it reproduces 
fairly well the limit of qmax cxd . The Bethe-Salpeter equation in hnite volume, can 
be written as, 

f-^ = V-^ -G (8) 


The energy levels in the box in the presence of interaction V correspond to the 
condition 

det{I - VG) = 0. (9) 

In a single channel, Eq. (|9]) leads to poles in the T amplitude when V~^ = G. 
Furthermore, for one channel, we can write the amplitude in inhnite volume T for 
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Figure 2: Left: Representation of 5G = G—G for D~^D*~ in function of qmax for ^/s = 3850 
MeV. The thick line represents the average of 5G for different values of qmax between 
1500 and 2500 MeV. Right: G'(solid) and (dashed) energy dependence of D^D*~ for 
LrriT^ = 2.0. Black dots correspond to energies {E = P^) where V~^ = G. Vertical dotted 
lines are the free energies in the box for DD*. 

the energy levels {Ei) as 

T = {G{E,) - G{Ei))-\ (10) 


4 Two channel case 

In the work [1] (see also DEI). a pole at ^/s = 3871.6 MeV is obtained using a 
subtraction constant of an = —1.265, with a binding energy of 0.2 MeV with respect 
to the neutral channel. When we address the inverse problem in the next section, for 
the sake of simplicity, we take only two channels, D'^D*~ and D^D , reevaluating the 
coupled channel calculation explained in section II (see Tabled]). Then, a new value 
an = —1.153 is needed in order to get the same position of the pole. The novelty 
of this study is the inclusion of two channel in the hnite box, where the energies are 
found using the condition of Eq. (jHI). As one can see in Fig. [3] we have two curves for 
each level, when for a single channel we had only a trajectory of the energy for each 
level. This feature is understood looking into Fig. [S] where the free energies for the 
channels D^D*~ and D^D are the dotted lines. Since the determinant of Eq. [9|has 
a zero between two asymptotes, the number of bound states in the box is doubled. 

5 The inverse problem 

Once we have determined the dependence of poles of T with L using the potential for 
the DD*, we want now to study the inverse problem. The idea is that QCD lattice 
data can be used to determine bound states of the DD system. For this purpose 
we assume that the lattice data are some discrete points on the energy trajectories 
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Figure 3: Left: L dependence of the energies of the poles for the two first levels of D~^D*~ 
and D^D* . Dotted lines correspond to the free energies. Right: Fit to the data. Dots with 
error bar are the synthetic data generated as explained in the text. Solid lines show the 
results obtained using the potential fitted to the synthetic data. 

obtained by us. Thus, we generate a set of data for some L for a value of the 
subtraction constant a = —1.153. In this case we generate 5 points in a range of 
LtUtt = [1.5, 3.5] and take 4 levels, this corresponds to n=0 and 1 in the momentum 
for both channels D~^D*~ and D^D . In addition, we simulate uncertainties in the 
obtained data, moving randomly by 1 MeV the centroid of the energies, then we 
assign an error of 2 MeV to these data. In Fig. [3] we show the simulated set of data. 
The second step is to choose the potential. We have chosen a potential with linear 
dependence in i/i. This is given by 

Vi = ai + bi ( 11 ) 

where \/s^ = m^io + m^o* is the energy of the hrst threshold, and i=l, 2 and 3 are the 

indices for each channel {i=l for D'^D*~ ^ i=2 for D^D and i=3 for the nondiagonal 
potential). Therefore, there are six parameters to determine in the potential. With all 
these ingredients, we do the £t, evaluating those values of the parameters in Eq. fill II 
that minimize the function. The error band is obtained in the standard method [7] 
varying randomly the parameters of the potential in a moderate range (10% change) 
and choosing the set of parameters that satisfy the condition < Xmm + 1- With 
these sets of parameters we determine the binding energy of the system with its 
dispersion from the pole of T = In both T and T we need a value of 

a to determinate G or G^^. The interesting thing that we observe is that the results 
for the binding energy are essentially independent of the choice of a. Changes in a 
revert on changes of V that compensate for it. We made choices of an between -1.2 
and -2.2. 


5 












B 

P 

Data 

AE 

AC 

ai 

0-2 

Parameters 

as 6i 


^3 


Results 

Pole Mean Pole 

(7 

4 

5 

2 

1 

-140.18 

-112.08 

-132.81 

-0.310 

0.074 

0.012 

2.32 

3871.51 

3871.49 

0.07 

4 

5 

5 

2 

-140.18 

-112.08 

-132.81 

-0.310 

0.074 

0.012 

0.79 

3871.51 

3871.25 

0.38 

4 

3 

2 

1 

-133.01 

-131.92 

-124.60 

-0.242 

0.048 

-0.075 

1.02 

3871.44 

3871.49 

0.18 

4 

3 

5 

2 

-120.09 

-98.19 

-150.94 

-0.377 

-0.075 

0.102 

0.28 

3871.41 

3871.15 

0.49 

2 

5 

2 

1 

-176.08 

-154.11 

-89.26 

9.92 

7.01 

-8.72 

0.259 

3871.70 

3871.47 

0.30 

2 

5 

5 

2 

-158.49 

-152.15 

-103.23 

4.56 

6.58 

-6.74 

0.982 

3871.34 

3871.30 

0.43 

2 

3 

2 

1 

-132.74 

-176.62 

-105.53 

3.23 

0.84 

-3.36 

0.074 

3870.51 

3870.48 

0.61 

2 

3 

5 

2 

-226.57 

-194.51 

-32.74 

31.81 

13.28 

-18.89 

0.942 

3869.49 

3870.37 

1.06 


Table 2: All possible set up changing number of branches {B), number of points (P), energy 
error bar (AE) and centroid of the energies (AC) and their set of parameters fitted. The 
columns denoted as Results are the obtained in the fit, the pole is determined with the 
parameters, and the mean pole and the dispersion are calculated as explained in the text. 
The results are for a = —1.25. As noted in the text, the use of different values of a change 
the potential but not the binding energy. Note that we quote values of total not the 
reduced one, which is always much smaller than 1. 

6 Results 

The results of the fits are shown in Table [2l where the hrst four columns determine 
the chosen set up of the synthetic data. The next columns are the fitted parameters, 
value of and pole position. The energy values in the “Pole” column correspond 
to the pole positions of the T matrix using the loop function evaluated with 
dimensional regularization together with the parametrized potential of Eq. flTT]) . To 
test the stability of the pole with the parameters, we vary randomly the parameters 
by 10%. If the new x^ calculated with those parameters is less than the x^ obtained 
in the fit plus one, we determine the pole position, otherwise it will be discarded. 
We iterate several times until we get 20 or 30 values of the pole positions. Then, we 
calculate the mean value of those pole positions and their dispersion a. The results 
are in the line with one should expect: fewer branches, fewer points or bigger errors 
which reverts into a higher dispersion in the binding energy. Since it is difficult for 
Lattice simulations to calculate higher levels, we have done also the test for the first 
level of energies for both channels, and in all cases the dispersion of the pole is higher 
than in the case where two levels are taken into account. Since the experimental 
errors in the binding of the X(3872) are of the order of 0.20 MeV, the exercise done 
is telling the level of precision demanded for the Lattice data if the experimental 
precision is to be matched. 


7 Conclusions 

We have studied the X(3872) state using coupled channels D^D*~ and D^D in a 
hnite box. This is done for a small binding energy. We have observed that in order 
to get a good precision in the binding energy, one does not need to extract the lattice 
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data with very small errors. Indeed, even with errors in the data points of 5 MeV, one 
can obtain the binding energy with 1 MeV (or even smaller value) precision. From a 
practical point of view, knowing that it is difficult to get four levels in actual Lattice 
calculations, it is rewarding to see that with only two levels one can get quite an 
accurate value for the binding, provided the levels are evaluated at several values of 
L with enough precision. We hope that this work gives a reference in the study of 
Lattice QCD for best strategies in order to obtain optimum values of the binding of 
the X(3872) state. 
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